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Abstract 

The human dental follicle partially differentiates into the periodontal ligament (PDL), but their biological functions are 
different. The gene-expression profiles of the dental follicle and PDL were compared using the cDNA microarray 
technique. Microarray analysis identified 490 genes with a twofold or greater difference in expression, 365 and 125 of 
which were more abundant in the dental follicle and PDL, respectively. The most strongly expressed genes in the 
dental follicle were those related to bone development and remodeling (EGFL6, MMP8, FRZB, and NELL1), 
apoptosis and chemotaxis (Nox4, CXCL13, and CCL2), and tooth and embryo development (WNT2, PAX3, FGF7, 
AMBN, AMTN, and SLC4A4), while in the PDL it was the tumor-suppressor gene WIF1. Genes related to bone 
development and remodeling (STMN2, IBSP, BMP8A, BGLAP, ACP5, OPN, BMP3, and TM7SF4) and wound 
healing (IL1, IL8, MMP3, and MMP9) were also more strongly expressed in the PDL than in the dental follicle. In 
selected genes, a comparison among cDNA microarray, real-time reverse-transcription polymerase chain reaction, 
and immunohistochemical staining confirmed similar relative gene expressions. The gene-expression profiles 
presented here identify candidate genes that may enable differentiation between the dental follicle and PDL. 
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Introduction 

The human dental follicle or dental sac is a loose 
ectomesenchymally derived connective tissue that surrounds 
the enamel organ and dental papilla of the developing tooth 
germ prior to eruption[1]. One of the biological functions of the 
dental follicle is the coordination of tooth eruption. Moreover, 
this tissue harbors progenitor cells for the periodontium[2], 
which, as the supporting tissue of the tooth, is composed of the 
periodontal ligament (PDL), alveolar bone, and the cementum 
that covers the tooth root surfaces[3]. The PDL tissue functions 
as a connector between the tooth and the alveolar jaw bone in 
the area surrounding the root surfaces[4]. Since the PDL 
originates from the dental follicle in tooth development, 
comparative gene-expression analyses may identify the unique 
features and specific functions of the two other dental tissues 
(i.e., alveolar bone and cementum). 



The introduction of cDNA microarray technology, which has 
become one of the standard laboratory tools for gene- 
expression analysis, has made it possible to analyze 
thousands of genes simultaneously. This technique has been 
used in the following three ways in dental research; 

1. To verify genes related to diseases or abnormal states, 
such as ameloblastoma^], dental caries[6,7], and the PDL 
under mechanical stress in vitro[8]. 

2. To compare two similar normal tissues in order to 
determine subtle genetic differences, such as permanent PDL/ 
deciduous PDL[8], fibroblast in the PDL/gingiva in v/fro[9]. 

3. To compare two different tissues that share the same 
origin, such as odontoblast/pulp tissue[10]. 

Such research has revealed the gene expression of the 
human tooth. However, no comparative study has yet been 
made of the dental follicle and PDL. 
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The present study evaluated and compared the gene- 
expression profiles of the dental follicle and PDL in humans 
with the aim of explaining their functions, such as eruption 
coordination and stress resorption. This information may be 
applicable to clinical problems such as eruption disturbance, 
and to periodontal tissue engineering. 

Materials and Methods 

1. Samples 

The experimental protocol was approved by the Institutional 
Review Board of Yonsei University Dental Hospital, and written 
informed consent to participate in the study was obtained from 
all of the subjects and their parents (#2-2011-0055). PDL 
samples were obtained from healthy permanent premolars 
(n=11; 1 male and 4 females, aged 11-19 years) that had been 
extracted for orthodontic reasons. Dental follicles were 
obtained from children (n=4; 3 males and 1 female, aged 6- 
7 years) during the extraction of supernumerary teeth. The 
extracted teeth and dental follicles were immediately frozen 
and stored in liquid nitrogen. The PDL tissues were obtained 
carefully using sterile curettes from the middle-third. The PDL 
tissues and dental follicles were then immediately submerged 
in RLT buffer, the component of the RNeasy Fibrous Mini kit 
(Qiagen, CA, USA). 

2. RNA isolation 

The dental follicle and PDL tissues were homogenized 
separately using a Bullet Blender Bead (Next Advanced, NY, 
USA), and total RNA was extracted from the homogenates 
using the RNeasy Fibrous Mini kit (Qiagen) in accordance with 
the manufacturer's instructions and as described previously[8]. 
The extracted RNA was eluted in 25 pi of sterile water and 
RNA concentrations were determined from absorbance values 
obtained at a wavelength of 260 nm as measured with the aid 
of a spectrophotometer (NanoDrop ND-1000, Thermo 
Scientific, IL, USA). The RNA samples used in this study had 
260/280 nm ratios equal to or greater than 1 .8. 

3. cDNA microarray analysis 

In the present study, global gene-expression analyses were 
performed using Affymetrix GeneChip Human Gene 1.0 ST 
oligonucleotide arrays (Affymetrix, CA, USA). The sample 
preparation was conducted in accordance with the 
manufacturer's instructions and recommendations. The 
average amount of RNA isolated from dental follicles or PDL 
tissues was 1 pg. RNA quality was assessed with the aid of an 
Agilent 2100 bioanalyzer using the RNA 6000 Nano Chip 
system (Agilent Technologies, Amstelveen, The Netherlands), 
and the quantity was determined using the NanoDrop ND-1000 
device (Thermo Scientific). 

As recommended by the manufacturer's protocol, 300 ng of 
total RNA from each sample was converted to double-strand 
cDNA. Using a random hexamer incorporating a T7 promoter, 
amplified RNA (cRNA) was generated from the double- 
stranded cDNA template though an in vitro transcription 
reaction and purified with the Affymetrix sample cleanup 



module. cDNA was regenerated through a random-primed 
reverse transcription (RT) using a dNTP mix containing dUTP. 
The cDNA was then fragmented by uracil-DNA glycosylase and 
apurinic/apyrimidinic endonuclease and restriction 
endonucleases, and end-labeled by terminal transferase 
reaction incorporating a biotinylated dideoxynucleotide. 
Fragmented end-labeled cDNA was hybridized to the 
GeneChip Human Gene 1.0 ST arrays for 16 hours at 45°C 
and 60 rpm, as described in the GeneChip Whole Transcript 
Sense Target Labeling Assay manual (Affymetrix). After 
hybridization, the chips were stained and washed in a 
GeneChip Fluidics Station 450 (Affymetrix) and scanned using 
a GeneChip Array scanner 3000 G7 (Affymetrix); the image 
data were extracted using Affymetrix Command Console 
software version 1.1 (Affymetrix). The raw file generated 
through this procedure contained expression intensity data and 
was used for the next step of the analysis. 

4. Microarray data analysis 

The expression data were generated by Affymetrix 
Expression Console software version 1.1 (Affymetrix). The 
Robust Multi-array Average (RMA) algorithm implemented in 
the Affymetrix Expression Console software was used to 
normalize the data. A one-way ANOVA was performed on the 
RMA expression values to determine whether genes were 
differentially expressed between the three groups. A multiple- 
testing correction was applied to the p values of the F-statistics 
to adjust the false discovery rate (Benjamini and Hochberg, 
1995). Genes with adjusted F-statistic p values of <0.05 were 
extracted. Microarray analysis identified 490 genes with 
expression differences of twofold or greater, 365 and 125 of 
which were more abundant in the dental follicle and PDL 
tissue, respectively. However, only strongly expressed genes in 
the dental follicles or PDL that differed by more than 4- or 2.5- 
fold from the signal of the control and each test group, 
respectively, were selected for further study. In order to classify 
the coexpression gene group with a similar expression pattern, 
hierarchical clustering and K-mean clustering were performed 
using MultiExperiment Viewer software version 4.4 
(www.tm4.org; Dana-Farber Cancer Institute, MA, USA). The 
Web-based tool Database for Annotation, Visualization, and 
Integrated Discovery (DAVID) was used for the biological 
interpretation of differentially expressed genes. These genes 
were classified based on data on gene function in the gene 
ontology of the Kyoto Encyclopedia of Genes and Genomes 
(KEGG) pathway database ( http://david.abcc.ncifcrf.gov/ 
home.jsp ). 

This microarray data set was approved by the Gene 
Expression Omnibus ( http://www.ncbi.nlm.gov/qeo/ ); its GEO 
accession number is GSE51342. 

5. Quantitative reverse-transcription polymerase chain 
reaction 

The single-stranded cDNA required for the polymerase chain 
reaction (PCR) analysis was produced using 500 ng of 
extracted total RNA as a template for RT (Superscript III 
Reverse Transcriptase and random primer, Invitrogen, UK). 
The RT reaction was performed at 65°C for 5 minutes, followed 
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by 25°C for 5 minutes, 50°C for 1 hour, and 70°C for 
15 minutes to inactivate the activity of the reverse 
transcriptase. The synthesized cDNA was diluted at 10:1 in 
distilled water and used as a template for quantitative RT-PCR, 
which was performed using the ABI7300 RT-PCR system 
(Applied Biosystems, Warrington, UK). Samples of 25 pi 
containing 1* Universal TaqMan Master Mix (4369016, Applied 
Biosystems), PCR primers at a concentration of 0.9 uM, and 
the diluted cDNA were prepared in triplicate. The amplification 
conditions were 50°C for 2 minutes and 95°C for 10 minutes 
followed by 40 cycles of 95°C for 15 seconds and 60°C for 
1 minute. The following TaqMan gene-expression assay 
primers (Applied Biosystems) were used: AMTN, CD36, 
CXCL13, DMP1, EGFL6, MMP8, MMP9, WIF1, and 18S rRNA. 
ABI 7300 SDS 1.3.1 software (Applied Biosystems) was used 
to record the fluorescence intensity of the reporter and 
quencher dyes; the results are plotted versus time, quantified 
as the cycle number. A precise quantification of the initial target 
was obtained by examining the amplification plots during the 
early log phase of product accumulation above background 
[the threshold cycle (Ct) number]. Ct values were subsequently 
used to determine ACt values (ACt=Ct of the gene minus Ct of 
the 18S rRNA control), and differences in Ct values were used 
to quantify the relative amount of PCR product, expressed as 
the relative change by applying the equation 2~ ACt . 

6. Immunohistochemical staining 

For immunohistochemical (IHC) staining, permanent teeth 
were fixed in 10% buffered formalin (Sigma, MO, USA) for 
1 day and then decalcified with 10% EDTA (pH 7.4; Fisher 
Scientific, TX, USA) for 8 weeks. The dental follicles and the 
permanent teeth were embedded in paraffin, and sectioned at 
a thickness of 3 pm. Specimens were subjected to IHC staining 
with antihuman AMTN (rabbit polyclonal, diluted 1:100; 
AM 2231 2, Abeam, Cambridge, UK), CXCL13 (BCA1, rabbit 
polyclonal, diluted 1:200; Ab1 12521, Abeam), DMP1 (rabbit 
polyclonal, diluted 1:100; Ab82351, Abeam), WIF1 (rabbit 
polyclonal, diluted 1:100; Ab71204, Abeam), and MMP9 (rabbit 
polyclonal, diluted 1:800; Ab38898, Abeam). The endogenous 
peroxidase activity was quenched by the addition of 3% 
hydrogen peroxide. Sections were incubated in 5% bovine 
serum albumin to block nonspecific binding. The primary 
antibodies in which the sections were incubated overnight were 
diluted so as to optimize the staining. After incubation, ready to 
use EnVision+ System HRP-labeled polymer antirabbit (K4003, 
ready to use; Dako North America, CA, USA) was applied for 
20 minutes. Color development was performed using labeled 
streptavidin biotin kits (Dako) following the manufacturer's 
instructions. The sections were counterstained with Gill's 
hematoxylin (Sigma). Control sections were treated in the 
same manner except for the treatment with primary antibodies. 

Results 

1. Dental follicle and PDL gene-expression profiles 

Complementary DNA microarray technology was used to 
compare multiple gene-expression profiles representative of 
the dental follicle and PDL tissues. The results indicated that 



the expression levels of 490 out of 33,297 (1 .49%) genes had 
changed by at least twofold in one tissue type relative to the 
other. In the dental follicle, the expression levels of 365 genes 
were at least double those in the PDL, while in the latter, the 
expression levels of 125 genes were at least double those in 
dental follicle tissue. The total data distribution and frequency 
were confirmed in density and box plots and in scatter plots of 
the standardized log intensity ratio versus the average intensity 
(Figure S1). Ultimately, 108 genes were analyzed further. In 
dental follicle tissue, 55 genes were up-regulated by at least 
fourfold relative to PDL tissue (Table 1 ), while 53 genes were 
up-regulated by at least 2.5-fold in the PDL relative to dental 
follicle tissue (Table 2). 

2. Gene-ontology analysis 

In order to identify the biological functions and features of the 
selected genes, the expressed data sets were organized into 
Gene Ontology Consortium (GO) grouping using the DAVID 
Web-based tool. These genes were then classified based on 
information regarding gene function using gene ontology from 
the KEGG pathway database. Figure 1 show GO classes with 
F-statistic p values of <0.05 for the two data sets analyzed. 

Biological process analysis revealed that genes related to 
cell surface receptor-linked signal transduction, cell adhesion, 
biological adhesion, neurologic system processes, and Wnt 
receptor signaling pathway were more strongly expressed in 
dental follicle tissue than in the PDL. Moreover, molecular 
function analysis revealed that genes related to ion binding, 
cation binding, metal ion binding, and calcium ion binding were 
more strongly expressed in dental follicle tissue than in the 
PDL. 

3. Quantitative RT-PCR 

Quantitative RT-PCR analysis was performed to verify the 
results obtained through cDNA microarray. The eight genes 
with expression levels differing by at least twofold between 
dental follicle tissue and the PDL were selected (Table S1). For 
the statistical analysis, the unpaired f-test (a<0.05) was 
performed to correlate the relative change with differential 
expression as detected by PCR. AMTN, EGFL-6, CXCL13, and 
MMP8 were up-regulated in dental follicle tissue relative to the 
PDL, while DMP1, WIF1 and CD36 were up-regulated in PDL 
tissue relative to the dental follicle (Table 3). The expression of 
MMP9 did not differ between dental follicle tissue and PDL 
tissue. These results are consistent with the microarray results. 

4. IHC staining 

The following five genes were the targets of the IHC study: 
AMTN, CXCL13, DMP1, WIF1, and MMP9 (Figure 2). AMTN 
staining was found only in the reduced enamel epithelium of 
the dental follicle. CXCL13 was broadly stained in the outer 
area of the dental follicles, and in Hertwig's epithelial root 
sheath (HERS). DMP1 was not stained in dental follicle tissue, 
but was found around the cementoblast layer. Similarly, WIF1 
was not stained in dental follicle tissue but was strongly stained 
in permanent PDL tissues, especially in the cementoblast layer. 
MMP9 was broadly stained in the outer area of the dental 
follicle but also in all layers of permanent PDL tissue. These 
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Table 1. Up-regulated genes in the dental follicle tissue (compared to periodontal ligament tissue). 



Name 


Gene symbol 


Fold change Gene Accession 


Cytoband 


ameloblastin (enamel matrix protein) 


AMBN 


94.88 


NM016519 


4q21 


amelotin 


AMTN 


22.23 


NM 2 12557 


4q13.3 


EGF-like-domain, multiple 6 


EGFL6 


17.76 


N M 0 15507 


Xp22 


microfibrillar associated protein 5 


MFAP5 


14.29 


NM 003480 


12p13.1-p12.3 


leucine-hch repeat-containing G protein-coupled receptor 5 


LGR5 


12.93 


NM 003667 


12q22-q23 


neuronal growth regulator 1 


NEGRI 


12.85 


NM_1 73808 


1p31.1 


coagulation factor II (thrombin) receptor-like 2 


F2RL2 


10.99 


NM_004101 


5q13 


hydroxyprostaglandin dehydrogenase 15-(NAD) 


HPGD 


10.66 


NM 000860 


4q34-q35 


gamma-aminobutyric acid (GABA) A receptor, alpha 3 


GABRA3 


9.60 


NM 000808 


Xq28 


PDZ domain containing ring finger 4 


PDZRN4 


8.61 


N M 0 13377 


12q12 


frizzled-related protein 


FRZB 


8.26 


NM 001463 


2qter 


statherin 


STATH 


7.57 


NM 003154 


4q13.3 


triadin 


TRDN 


7.31 


NM 006073 


6q22.31 


chemokine (C-X-C motif) ligand 13 


CXCL13 


7.21 


NM 006419 


4q21 


leucine rich repeat neuronal 1 


LRRN1 


7.10 


NM 020873 


3p26.2 


C1q and tumor necrosis factor related protein 3 


C1QTNF3 


7.09 


NMJ81435 


5p13 


spondin 1 , extracellular matrix protein 


SPON1 


6.66 


NM 006108 


11p15.2 


ectonucleotidepyrophosphatase/phosphodiesterase 5 (putative) 


ENPP5 


6.52 


NM 021572 


6p21.1-p11.2 


epithelial cell adhesion molecule 


EPCAM 


6.41 


NM 002354 


2p21 


solute carrier family 4, sodium bicarbonate cotransporter, member 4 


SLC4A4 


6.33 


NM 001 098484 


4q21 


ALX homeobox 1 


ALX1 


6.29 


NM 006982 


12q21.31 


NADPH oxidase 4 


NOX4 


6.23 


N M 0 16931 


11q14.2-q21 


fibroblast growth factor 7 


FGF7 


5.75 


NM 002009 


15q21.2 


matrix metallopeptidase 8 (neutrophil collagenase) 


MMP8 


5.58 


NM 002424 


11q22.3 


small nucleolar RNA, C/D box 114-26 


SNORD 114-26 


5.58 


NR003219 


14q32 


sarcoglycan, zeta 


SGCZ 


5.49 


NMJ39167 


8p22 


5.38 


hyaluronan and proteoglycan link protein 1 


HAPLN1 


5.12 


NM 001884 


5q14.3 


CD24 molecule 


CD24 


4.94 


N M 0 13230 


6q21 


transient receptor potential cation channel, subfamily A, member 1 


TRPA1 


4.86 


NM 007332 


8q13 


carboxypeptidase X (M14 family), member 2 


CPXM2 


4.85 


NMJ98148 


10q26.13 


paternally expressed 10 


PEG10 


4.75 


N M 0 15068 


7q21 


sema domain, immunoglobulin domain (Ig), short basic domain, secreted, (semaphorin) 3A 


SEMA3A 


4.67 


NM 006080 


7p12.1 


flavin containing monooxygenase 3 


FM03 


4.66 


NM 006894 


1q24.3 


wingless-type MMTV integration site family member 2 


WNT2 


4.65 


NM_003391 


7q31.2 


thrombospondin 4 


THBS4 


4.62 


NM 003248 


5q13 


pappalysin 2 


PAPPA2 


4.59 


NM 020318 


1q23-q25 


G protein-coupled receptor 64 


GPR64 


4.53 


NM 001 079858 


Xp22.13 


elongation of very long chain fatty acids (FEN1/Elo2, SUR4/Elo3, yeast)-like 2 


ELOVL2 


4.51 


N M 0 17770 


6p24.2 


calcium channel, voltage-dependent, alpha 2/delta subunit 3 


CACNA2D3 


4.48 


NM 018398 


3p21.1 


complement component 3 


C3 


4.46 


NM 000064 


19p13.3-p13.2 


collagen, type XIV, alpha 1 


COL14A1 


4.36 


NM021110 


8q23 


Ras protein-specific guanine nucleotide-releasing factor 2 


RASGRF2 


4.31 


NM 006909 


5q13 


activated leukocyte cell adhesion molecule 


ALCAM 


4.30 


NM 001627 


3q13.1 


neural cell adhesion molecule 2 


NCAM2 


4.23 


NM 004540 


21q21.1 


dickkopf homolog 2 (Xenopuslaevis) 


DKK2 


4.22 


N M 0 14421 


4q25 


hemicentin 1 


HMCN1 


4.20 


NM 031935 


1q25.3-q31.1 


core 1 synthase, glycoprotein-N-acetylgalactosamine 3-beta-galactosyltransferase, 1 


C1GALT1 


4.17 


NM 020156 


7p14-p13 


NEL-like 1 (chicken) 


NELL1 


4.16 


NM 006157 


11p15.1 


chromosome 12 open reading frame 75 


C12orf75 


4.14 


NM 001 145199 


12q23.3 


odz, odd Oz/ten-m homolog 2 (Drosophila) 


ODZ2 


4.13 


AK302302 


5q34 


paired box 3 


PAX3 


4.09 


NMJ81458 


2q35-q37|2q35 


solute carrier family 38, member 1 


SLC38A1 


4.08 


NM_030674 


12q13.11 


gremlin 1 


GREM1 


4.04 


N M 0 13372 


15q13.3 
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Table 1 (continued). 





Name 


Gene symbol 


Fold change Gene Accession 


Cytoband 


solute carrier family 1 (glial high affinity glutamate transporter), member 3 


SLC1A3 


4.01 


NM 004172 


5p13 


Kallmann syndrome 1 sequence 


KAL1 


4.00 


NM 000216 


Xp22.32 


microfibrillar-associated protein 4 


MFAP4 


3.97 


NM_001 198695 


17p11.2 


dermatan sulfate epimerase-like 


DSEL 


3.96 


NM 032160 


18q22.1 


laminin, alpha 2 


LAMA2 


3.95 


NM 000426 


6q22-q23 


vesicle amine transport protein 1 homolog (T. californica)-like 


VAT1L 


3.91 


NM 020927 


16q23.1 


chemokine (C-C motif) ligand 2 


CCL2 


3.88 


NM_002982 


17q11.2-q12 



doi: 10.1371/journal.pone.0084201.t001 



Wnt receptor signaling. 

cell death 
blood vessel development 
cell motion 
neurologic system process 
cell-cell signaling 
transcription 
response to wound 
cell-surface- receptor-linked, 
biological adhesion 
cell adhesion 

a. 




nucleoside binding 
cell morphogenesis 
DNA binding 
growth factor binding 
ion channel activity 
cytokine activity 
fibroblast growth, 
structural molecule.. 
calcoum ion binding 
metal-ion binding 
cation binding 




■ periodontal 
ligament 



■ dental 
follicle 



Figure 1. Gene-ontology analysis of the periodontal ligament(PDL) and the dental follicle, a. Main categories of genes 
expressed specifically in the dental follicle and periodontal ligament (PDL) tissues on the basis of their biological processes, b. Main 
categories of genes expressed specifically in dental follicle and PDL tissues on the basis of their molecular function (x-axis: number 
of involved genes, F-statistic p<0.05). 

doi: 10.1371/journal.pone.0084201.g001 



results are consistent with those of the cDNA microarray 
analysis at the protein level. 

Discussion 

The cDNA microarray results indicate that the expression 
levels of 490 out of 33,297 (1.49%) genes had changed by at 
least twofold in one tissue type relative to the other. Although 
this finding is not dissimilar to that of McLachlan et al.[7] in 
dental caries (445/15,000, 2.96%), Song [8] reported that only 
51 of >20,000 genes (<0.25%) were differentially expressed 
between the PDLs of permanent teeth and deciduous teeth. 
Differences in the homogeneity between dental follicle and PDL 
tissue from different tooth types may explain this discrepancy, 
with PDL tissues from different tooth types possibly being more 
homogeneous than the dental follicle. 

In general, genes associated with bone development and 
remodeling (EGFL6, MMP8, FRZB, and NELL1), apoptosis and 
chemotaxis (Nox4, CXCL13, and CCL2), and tooth and embryo 
development (WNT2, PAX3, FGF7, AMBN, AMTN, and 
SLC4A4) were more strongly expressed in the dental follicle 
than in the PDL. These genes are thought to be associated 
with the biological function of dental follicle tissue, such as 
eruption coordination. 



In detail, bone remodeling requires an intimate cross-talk 
between osteoclasts and osteoblasts. EGFL6, a member of the 
epidermal growth factor repeat superfamily of proteins, induces 
angiogenesis via a paracrine mechanism in which EGFL6 is 
expressed in osteoblast-like cells but promotes the migration 
and angiogenesis of endothelial cells [1 1], MMP8 is involved in 
the breakdown of extracellular matrix in normal physiological 
processes, such as tissue remodeling[12], FRZB (frizzled motif 
associated with bone development) plays an important role in 
bone development [13]. NELL1 promotes osteoblast cell 
differentiation and mineralization. Runx2, Osterix, and NELL1 
promote crosstalk during osteogenesis [14], 

Nox4 promotes the apoptosis of oral epithelium at the 
eruption site and so is thought to be related to the tooth- 
eruption mechanism [15]. CXCL13 [chemokine (C-X-C motif) 
ligand 13] promotes the migration of B lymphocytes and is 
strongly expressed in the liver, spleen, and lymph nodes [16]. 
The number of CXCL13 cells is reportedly significantly higher 
in periodontitis, suggesting that this condition is associated with 
alveolar bone resorption. CCL2 (MCP-1) is a well-known 
chemokine for monocytes and a prime candidate for recruiting 
the osteoclast precursors into dental follicle tissue [17]. 

Genes associated with tooth and embryo development were 
also significantly more strongly expressed in the dental follicle 
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Table 2. Up-regulated genes in the periodontal ligament 
tissue (compared to dental follicle). 





Gene 


Fold 








Name 


symbol 


change 


Gene access 


cytoband 


WNT inhibitory factor 1 


WIF1 


12.63 


NM 


.007191 


12q14.3 


dentin matrix acidic 
phosphoprotein 1 


DMP1 


1 0 67 


NM 


nn4An7 

UU^f'fU / 


4q21 


acid phosphatase 5, 
tartrate resistant 


ACP5 


6.78 


NM 


001111035 


19p13.3- 
p13.2 


small proline-rich protein 
2A 


SPRR2A 


6.60 


NM 




1q21-q22 


bone morphogenetic 
protein 8a 


DlvlroA 


6 15 


NM 


A Q A onn 

I o i ouy 


1 p34 3 


transmembrane 4 L six 
family member 19 


TM4QF1Q 


6.12 


NM 


A QO/1 C A 

1 00401 


3q29 


tenascin N 


TNN 


6.02 


NM 


022093 


1q23-q24 


stathmin-like 2 


STMN2 


5.57 


NM 


007029 


8q21.13 


bone gamma- 












carboxyglutamate (gla) 


BGLAP 


5.10 


NM 


199173 


1q25-q31 


protein 












signal peptide, CUB 
domain, EGF-like 3 


SCUBE3 


4.55 


NM 


1 52753 


6p21.3 


CD36 molecule 












(thrombospondin 


CD36 


4.54 


NM 


001001548 


7q11.2 


receptor) 












ATPase, H+ transporting, 












lysosomal 38kDa, VO 


ATP6V0D2 


4.25 


NM 


152565 


... 


subunit d2 












small Cajal body-specific 
RNA17 


QPARNA17 


4 21 


NR 


003003 


18q21 1 


matrix metallopeptidase 9 












(gelatinase B, 92kDa 
gelatinase, 92kDa type IV 


MMP9 


4.20 


NM 


004994 


20q11.2- 
q13.1 


collagenase) 












multiple EGF-like- 
domains 10 


MEGF10 


4.14 


NM 


032446 


5q33 


interleukin 1, beta 


IL1B 


4.11 


NM. 


000576 


2q14 


G protein-coupled 
receptor 109B 


GPR109B 


4.10 


NM 


006018 


12q24.31 


interleukin 8 


IL8 


4.07 


NM 


000584 


4q13-q21 


3.98 


anoctamin 5 


AN05 


3.94 


NM 


213599 


11p14.3 


integrin-binding 
sialoprotein 


IBSP 


3.91 


NM 


004967 


4q21.1 


matrix metallopeptidase 3 












(stromelysin 1 , 


MMP3 


3.80 


NM 


002422 


11q22.3 


progelatinase) 












matrix extracellular 
phosphoglycoprotein 


MEPE 


3.73 


NM 


001184694 


4q21.1 


interferon, gamma- 
inducible protein 30 


IFI30 


3.61 


NM 


006332 


19p13.1 


creatine kinase, brain 


CKB 


3.46 


NM. 


.001823 


14q32 


integrin, alpha 10 


ITGA10 


3.42 


NM 


003637 


1q21 


small Cajal body-specific 
RNA 17 


SCARNA17 


3.38 


NR. 


003003 


18q21.1 



Table 2 (continued). 





Gene 


Fold 








Name 


symbol 


change 


Gene access 


cytoband 


BMP and activin 












membrane-bound 
inhibitor homolog 


BAMBI 


3.19 


NM 


Ul Z J4z 


10p12.3- 
p11.2 


(Xenopuslaevis) 












pannexin 3 


PANX3 


3.18 


NM. 


052959 


11q24.2 


solute carrier family 37 












(glycerol-3-phosphate 


SLC37A2 


3.18 


NM 


198277 


11q24.2 


transporter), member 2 












transcription factor AP-2 












beta (activating enhancer 


TFAP2B 


3.14 


NM 


.003221 


6p12 


binding protein 2 beta) 












bone morphogenetic 
protein 3 


BMP3 


3 09 


NM 


UU I Z.KJ I 


4q21 


distal-less homeobox 1 


DLX1 


2.98 


NM 


178120 


2q32 


interferon induced 
transmembrane protein 5 


IFITM5 


2.97 


NM 


001025295 


11p15.5 


bone morphogenetic 
protein 8b 


BMP8B 


2.94 


NM 


UU1 / Z{J 


1 p35-p32 


myelin protein zero 


MPZ 


2.93 


NM 


000530 


1q23.3 


MAM domain containing 
2 


MAMDC2 


2.91 


NM 


1 53267 


9q21.12 


dachsous 2 (Drosophila) 


DCHS2 


2.90 


NM 


017639 


4q31.3 


chondroadherin 


CHAD 


2.85 


NM 


001267 


17q21.33 


serpin peptidase inhibitor, 












clade E (nexin, 
plasminogen activator 


SERPINE1 


2.83 


NM 


000602 


7q21 .3- 


inhibitor type 1), member 
1 










selectin E 


SELE 


2.77 


NM 


000450 


1q22-q25 


transmembrane 7 
superfamily member 4 


TM7 C 1F4 
i i vi / or^t 


2 77 


NM 


UoU / oo 


8q23 


aquaporin 9 


AQP9 


2.76 


NM 


020980 


15q 


protein tyrosine 












phosphatase, receptor 












type, f polypeptide 


PPFIA2 


2.76 


NM 


003625 


12q21.31 


(PTPRF), interacting 












protein (liprin), alpha 2 












phospholipase A2, group 












IIA (platelets, synovial 


PLA2G2A 


2.75 


NM 


000300 


1p35 


fluid) 












lysophosphatidic acid 
receptor 3 


LPAR3 


2.71 


NM 


012152 


1 p22 3 


rhomboid, veinlet-like 2 
(Drosophila) 


RHBDL2 


2.69 


NM 


017821 


1p34.3 






2.67 








discs, large (Drosophila) 












homolog-associated 


DLGAP5 


2.64 


NM 


014750 


14q22.3 


protein 5 












major histocompatibility 












complex, class II, DQ 


HLA-DQA1 


2.60 


NM 


002122 


6p21.3 


alpha 1 












von Willebrand factor D 
and EGF domains 


VWDE 


2.58 


NM. 


.001135924 


7p21.3 



PLOS ONE | www.plosone.org 



6 



December 2013 | Volume 8 | Issue 12 | e84201 



Genetic Comparison of the Dental Follcle and PDL 



Table 2 (continued). 







Gene 


Fold 








Name 


symbol 


change 


Gene access 


cytoband 


family with sequence 
similarity 40, member B 


FAM40B 


2.57 


NM 


020704 


7q32.1 


chemokine (C-X-C motif) 
receptor 1 


CXCR1 


2.57 


NM 


000634 


2q35 


retinol binding protein 4, 
Dlasma 


RBP4 


2.52 


NM 


006744 


10q23-q24 


G protein-coupled 
receptor 109A 


GPR109A 


2.50 


NM 


177551 


12q24.31 


microRNA 487a 


MIR487A 


2.50 


NR 


030162 


14q32.31 



doi: 10.1371/journal.pone.0084201.t002 



Table 3. The relative difference in gene mRNA expression 
in dental follicle and periodontal ligament tissue. 





Gene 


Relative Gene Expression (mean±SD) 




Dental Follicle Tissues Periodontal Ligament Tissues 


amtn" 


2741.88±245.81 1 


EGFL-6" 


265.40±79.19 1 


CXCL13" 


12.15±1 .23 1 


MMP8" 


3.48±0.70 1 


DMP1" 


1 51.59±14.64 


WIF1" 


1 233.17±24.40 


CD36" 


1 5.32±1.26 


MMP9 


1 1.34±0.26 



"p<0.01 

doi: 10.1371/journal.pone.0084201.t003 



than in the PDL. WNT2 is a member of the WNT gene family, 
and the Wnt pathway is crucial for tooth development, 
embryogenesis, and odontoblast and ameloblast differentiation 
[18]. Previous studies found that Wnt signaling inhibits cell 
differentiation and promotes cell proliferation in cementoblasts 
in vitro [19]. Furthermore, the Wnt signaling pathways may 
participate in the formation of the tooth eruption passage [15]. 

The tumor-suppressor gene WIF1 was the most strongly 
expressed in the PDL. In addition, genes associated with 
wound healing (IL1, IL8, MMP3, and MMP9) and bone 
development and remodeling (STMN2, IBSP, BMP8A, BGLAP, 
ACP5, OPN, BMP3, TM7SF4) were up-regulated relative to the 
dental follicle. 

WIF1 was the most strongly expressed gene in the PDL, 
which means that, conversely, the Wnt signal would have been 
strongly expressed in dental follicle tissue. WNT2 and DKK2 
are associated with WNT signaling, and the genes encoding 
both were strongly expressed in the dental follicle in this study. 
The Wnt signaling pathway has been implicated in a wide 
range of biological process, from maintaining the 
pluripotentiality of stem cells to the induction of specific tissues 
and organs during development. WIF1 functions as a tumor- 
suppressor gene, and it has been shown to be epigenetically 
silenced in various cancers[20]. Genes related to wound 



healing were significantly strongly expressed. IL1 and IL8 are 
important mediators of the inflammatory response and play an 
important role in the induction of osteoclastic bone resorption 
during orthodontic tooth movement [21]. MMP3 and MMP9 are 
regulated by fibroblasts in the PDL and are induced in 
fibroblasts and immune cells that are recruited in the pathologic 
situation [22]. 

Genes associated with bone development and remodeling 
were significantly more strongly expressed in the PDL than in 
the dental follicle. STMN2 is a novel marker of osteogenesis 
and osteoblasts [23]. IBSP (integrin-binding sialoprotein) is a 
major structural protein of the bone matrix [24]. BMP8A is 
related to osteogenesis and skeletal development [24,25]. 
BGLAP (osteocalcin) was first identified in the mineralized 
matrix of bone and is associated with bone development [26]. 
ACP5 (TRAP) exhibits optimal activity in acidic conditions and 
is strongly expressed by osteoclasts, activated macrophages, 
and neurons. OPN and bone sialoprotein are strongly 
associated with TRAP when phosphorylated [27]. BMP3 is a 
major structural protein of the bone matrix and osteoclast 
activity in fractured sites. This gene is a member of the 
transforming growth factor beta superfamily and is localized to 
the latest stage of periodontium formation by inhibiting the 
signaling provided by other bone morphogenetic proteins [28]. 
TM7SF4 is related to osteoclastogenesis and the regulation of 
osteoclast fusion [29]. 

Quantitative RT-PCR analyses were carried out to verify our 
cDNA microarray results. AMTN, EGFL-6, CXCL13, and MMP8 
were up-regulated in dental follicle tissue relative to the PDL, 
while DMP1, WIF1 and CD36 were up-regulated in the PDL 
tissue relative to the dental follicle. These findings are 
consistent with the microarray results. 

To better understand the roles of the differentially expressed 
genes identified by the microarray analyses between dental 
follicle and PDL tissues in humans, their cellular origins were 
identified using IHC analysis. AMTN was found only in the 
reduced enamel epithelium of dental follicle tissue. CXCL13 
was broadly stained in the outer area of the dental follicle and 
in HERS. DMP1 was not stained in the dental follicle, but was 
found around the cementoblast layer. WIF1 was not stained in 
dental follicles, but was strongly stained in permanent PDL 
tissue, especially the cementoblast layer. MMP9 was broadly 
stained in the outer area of the dental follicles, but also strongly 
in all layers of permanent PDL tissue (Figure 2). These results 
are consistent with those of the cDNA microarray analysis at 
the protein level. 

In conclusion, genes associated with bone development and 
remodeling were strongly expressed in both PDL and dental 
follicle tissues, and in vivo would induce tooth eruption or tooth 
movement. The two tissues also expressed specific genes 
related to their functions. Interestingly, WIF1, which was 
strongly expressed in the PDL, is thought to represent a 
counter example of the pluripotentiality of the stem cells of the 
dental follicle. The findings of this study are expected to 
contribute to future research in the field of clinical problems 
such as eruption disturbance, and to research in periodontal 
tissue engineering by demonstrating the differences between 
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Figure 1 
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Figure 2. Immunohistochemical (IHC) staining of dental follicle and permanent PDL tissues. (A, F) IHC staining for AMTN in 
the dental follicle and (K, P) permanent PDL tissues. (B, G) IHC staining for CXCL13 in the dental follicle and (L, Q) permanent PDL 
tissue. (C, H) IHC staining for DMP1 in the dental follicle and (M, R) permanent PDL tissue. (D, I) IHC staining for WIF1 in the dental 
follicle and (N, S) permanent PDL tissue. (E, J) IHC staining for MMP9 in the dental follicle and (O, T) permanent PDL tissue. (Scale 
bars: 20 urn in A-E and P-T; 100 urn in F-O.) 

doi: 10.1371/journal.pone.0084201.g002 



the PDL and dental follicle tissues at the molecular biological 
level. 



Supporting Information 

Figure S1. Scatter plots of the differences in the log 
intensities versus the average log intensity between three 
dental follicle tissue samples and three PDL tissue 
samples. In each plot the x-axis (A) is 0.5 x [log 2 (case)+log 2 
(control)] and the y-axis (M) is log 2 (case/control). The data of 
all plots were normally distributed. 
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